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explored by means of a number of soluble models. It is found that in all cases the 
energy remains real for small values of the coupling constant, but becomes complex if 
the coupling becomes stronger than some critical value. For a quadratic non-Hermitian 
Ph, "PT-symmctric Hamiltonian coupled to an arbitrary real Hermitian "PT-symmetric 

pH ' Hamiltonian, the reality of the ground-state energy for small enough coupling constant 

is established up to second order in perturbation theory. 
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1. Introduction 



^ Since the work by Bender and Boettcher jlj on non-Hermitian but PT-symmetric 

Hamiltonians, subsequent research has gone through various stages. First came an 
exploration of various non-Hermitian generalizations of well-known soluble models to 
determine if their spectra is real. However, reality of the spectrum does not by itself 
guarantee a viable quantum theory. One also needs a probabilistic interpretation, and 
since the most obvious choice of metric for a PT-symmetric model is not positive 
definite, a Hilbert space endowed with this metric does not represent a physical 
framework for quantum mechanics. Instead, one must find an alternative, positive- 
definite metric [2j[3j|l], which is dynamically determined by the particular Hamiltonian 
in question. It was further shown [5] that this metric 77 = e~® provides a similarity 
transformation from the non-Hermitian H to an equivalent Hermitian H. This 
equivalent Hermitian Hamiltonian was subsequently constructed, often in perturbation 
theory only, in a variety of models P, [7J E] • 
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These investigations were all concerned with isolated non-Hermitian systems, but 
more recently much attention has been given to situations where a non-Hermitian 
system interacts with the world of Hermitian quantum mechanics. For example, Ref. [9] 
examined a non-Hermitian analogue of the Stern-Gerlach experiment in which the role 
of the intermediate inhomogeneous magnetic field flipping the spin is replaced by an 
apparatus described by a non-Hermitian Hamiltonian. This type of set-up has been 
further elaborated by Assis and Fring [10] and Gunther et al. [11], and it continues to 
be the focus of lively discussion [12], [131 El EE EH] . Recently, Ref. p2] explored the 
problem of scattering from localized non-Hermitian potentials. 

It is in this spirit that we investigate the nature of the energy spectrum when 
Hermitian and non-Hermitian systems, individually having real spectra, are coupled 
together. In Sec. [2] we first look at a simple matrix model and then in Sec. [3] we 
explore couplings of various non-Hermitian quadratic Hamiltonians to a simple harmonic 
oscillator. Section H] examines the reality of the spectrum in perturbation theory for a 
complex quadratic PT-symmetric Hamiltonian coupled to a generic real PT-symmetric 
and Hermitian Hamiltonian. We summarize our results in Sec. [51 



2. Simple Matrix Model 

In this section we consider coupling a Hermitian matrix Hamiltonian 



Hi = al + bai 



a b 
b a 



(a, b real) 



(1) 



to the non-Hermitian "PT-symmetric matrix Hamiltonian introduced in Ref. [2] 



H 2 = r(l cos 9 + i sm.60-3) + sai 



re 



re 



(r, s, 9 real), (2) 



where 1 is the identity matrix and o~k are the Pauli matrices. We choose the parameters 
r, s, and 8 so that the inequality s 2 > r 2 sin 2 6 is satisfied; this inequality guarantees 
that the eigenvalues of if 2 are real [2]. The parity operator is taken as 



V = 01 



1 

1 



(3) 



while T implements complex conjugation. Note that both Hi and H 2 are symmetric 
under VT, and each separately has real eigenvalues. 

To couple these two systems together we take the direct sum, but introduce nonzero 
elements in the off-diagonal sector: 
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These are chosen in such a way that H remains invariant under combined parity 
reflection and time reversal, where the parity operator is given by 



V = 1 ® (Ji 
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o\ 
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1° 
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°) 



(5) 



and time reversal is complex conjugation. 

The question is whether the eigenvalues of this combined system remain real, and if 
there is any constraint on the strength of the coupling parameter e. As a specific example 
we choose a = b = 1 in H±, r = s = 1 and 6 = tt/6 in H2, and A = C = 1, B = D = 
in the coupling matrices. Numerically we find that in this case the eigenvalues remain 
real until e exceeds a critical value of approximately 0.7045. For other choices of the 
parameter the picture is similar: In all cases the eigenvalues remain real for e sufficiently 
small. In some cases they appear to remain real for all values of e; in others, as in the 
example above, they first become complex but then become real again at a second 
critical value of e. 

3. Coupled Quadratic Hamiltonians 

To determine if the energy levels of a coupled system of Hamiltonians are real, our 
strategy here will be to find the Q operator by using the condition that 

#t = e -Q He Q , (6) 

and then to construct the equivalent Hermitian Hamiltonian H by using the similarity 
transformation 

H = e- Q / 2 He Q ' 2 . (7) 

In some cases the resulting H will need to be diagonalized by a further unitary 
transformation in order to identify the spectrum. 

3.1. Simple Harmonic Oscillator Coupled to a Shifted Simple Harmonic Oscillator 

In this subsection we consider a quantum system described by the interaction of a 
conventional and a "PT-symmetric Hamiltonian: 

H = (p 2 + x 2 ) + (q 2 + y 2 + 2iy) + 2exy. (8) 

Since the Hamiltonian ([8]) is quadratic, we expect the Q operator to be linear in the 
momentum variables, 

Q = ap + pq, (9) 

which will produce the coordinate shifts 

x — ► x + ia and y — > y + i(3. (10) 



Interactions of Hermitian and non-Hermitian Hamiltonians 



4 



We determine a and f3 by the condition which gives 

a = - 6 and (3 = - 9 . (11) 
1 — 1 — e z 

It is somewhat surprising that we are able to determine a and (3 because the condition 

jji _ e -Qff e Q translates into a system of three coupled linear equations with only two 

unknowns a and (3. Yet, there is a unique solution. However, note that the solution 

becomes singular as \e\ reaches 1. 

Given Q we construct H according to Eq. ([7]), which in this case produces 

H = e- Q/2 He Q/2 = p 2 + x 2 + q 2 + y 2 + 2exy + - 7 . (12) 

1 — e z 

We have identified the equivalent Hermitian Hamiltonian H, but it must still be 
diagonalized. To do so we change variables from p, q, x, and y to P, Q, X, and Y: 

p = aP + bQ, 
q = cP + dQ, 
x = eX + fY, 

y = gX + hY. (13) 

We determine the unknown coefficients a through h by requiring that (i) the canonical 
commutation relations 

[p, q] = 0, [x, y] = 0, [y,p] = 0, [x, q) = 0, [x,p] = i, [y, q] = i (14) 

be maintained, and that (ii) H, when expressed in terms of X, Y, P, and Q contains 
no crossterms. 

These two sets of conditions translate into six equations for the coefficients. The 
solutions to these six equations are: 

1 1 T] T] 

c = ria, d = —rib, e = — , f = — , q = — , h = — -, (15) 
' ' ' ' 2a J 26' y 2a 2b' V ; 

where rj = ±1 and a and b are arbitrary. The resulting H is given by 

H = 2a 2 P 2 + l -^X 2 + 2b 2 Q 2 + ^r^Y 2 + — l —, (16) 
2cr 2b z 1 — e z 

which is the sum of two uncoupled quantum-harmonic-oscillator Hamiltonians. 

Since the energy levels of the general quantum harmonic oscillator Hamiltonian 

H = Ap 2 + Bx 2 are 

E n = (2n + 1)VAB (n = 0,1,2,...), (17) 
the energy levels of H in fflBT) are 

E m>n = (2m+l)VTT~e + (2n+l)VT^l+ (18) 

This result is independent of the constants a and b as well as the choice of sign of rj. The 
energy diverges at the critical value |e| = 1 and becomes complex for |e| > 1. Thus, there 
are two regions, depending on whether the Hermitian component of the Hamiltonian is 
coupled strongly or weakly to the non-Hermitian component of the Hamiltonian. 
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This result is not specific to the choice of coefficients in (jS]), as we now show by 
considering the more general Hamiltonian 

H = {p 2 + u\x 2 ) + {q 2 + u 2 2 y 2 + 2i\y) + 2exy. (19) 

Again, we take the Q operator to have the form Q = ap + flq, and determine a and (5 
by the condition that H> = e~ Q He Q . We obtain 

a = l ,2 i t f 2 aIld P = - n 2 U ; 2 X f 2 - ( 20 ) 

The equivalent Hermitian Hamiltonian is then 

\ 2 u 2 

H = p 2 + u 2 x 2 + q 2 + u 2 y 2 + 2exy + 2 2 1 (21) 

Making the ansatz in ffTB"]) . we now obtain the unknown coefficients a through h by 
following the same procedure as above. The solutions are: 

A b 1 

c = 7a, a = , e 



7' a(l + 7 2 )' 

/= 6(1 +7 2 )' ^ = a(l + 7 2 )' /l = ~6(l+7 2 )' ^ 
where 7 satisfies the quadratic equation 

e7 2 + 2£> 7 - e = (23) 

and we define 

D=l((4-u%) and 5 = i(cj 2 + cu 2 ). (24) 
The resulting Hermitian Hamiltonian H is 

# = p 2 + fi 2 x 2 + Q 2 + fi 2 F 2 + ff 1 - , (25) 

where we have used the freedom in the choice of a and 6 to set a 2 = 1/(1 + 7 2 ) and 
b 2 = 7 2 / (1 + 7 2 ) and where the parameters Q\ and ^2 are given by 

Ql = S± VD 2 + e 2 and fi 2 = St VD 2 + e 2 . (26) 

The energy levels of the Hamiltonian (1251) are 



£ m , n = (2m + l)y/S + (D 2 + e 2 ) 1/2 

+ (2n + l)^S-(D 2 + e 2 f 2 + - . (27) 

Again, we find that the energy diverges, this time at the critical value |e| = 001002, and 
for |e| larger than this value the energy becomes complex. Thus, again there are two 
regions, depending on whether the Hermitian component of the Hamiltonian is coupled 
strongly or weakly to the non-Hermitian component of the Hamiltonian. 
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3.2. Two Coupled Shifted Simple Harmonic Oscillators 

The pattern that we observed in the previous subsection re-emerges when we consider 
two coupled PT-symmetric non-Hermitian Hamiltonians: 

H = {p 2 + x 2 + 2i\x) + {q 2 + y 2 + 2i/xy) + 2exy. (28) 



As before, we choose Q = ap + /3q, and determine a and (3 by condition ([6]). This gives 

eu — X , _ eX — u , n . 

a = ■£ — - and (3 = - £ . (29) 

1 — e z 1 — e z 

Applying (j7]), we obtain the equivalent Hermitian Hamiltonian 

H = e- Q ' 2 He Q ' 2 = p 2 + x 2 + q 2 + y 2 + 2exy + - + f- ~ ^ , (30) 

1 — e z 

which is exactly the same as the result in (1T2"1) . apart from the additive constant. 

3.3. Simple Harmonic Oscillator Coupled to Swanson Hamiltonian 

Here we consider the non-Hermitian system described by the Swanson Hamiltonian 
[IB] , written in terms of coordinate and momentum variables instead of creation and 
annihilation operators: 

H = {p 2 + x 2 ) + {q 2 + y 2 + ic{q, y} + ) + 2exy. (31) 

We can exploit the ambiguity in Q for the Swanson Hamiltonian itself to choose 
Q = —cy 2 . This shifts q — >• q — icy but leaves y, and hence the coupling term 2exy, 
unchanged. The equivalent Hermitian Hamiltonian is then 

H = p 2 + x 2 + q 2 + (1 + c 2 )y 2 + 2exy, (32) 

which can be diagonalized to give 

h = P 2 + ttjX 2 + Q 2 + £l 2 2 Y 2 , (33) 

where 



tt 2 lj2 = 1 + \c 2 [1 ± v /l + 4e 2 /c 4 J • (34) 

Notice that the eigenvalues 

E m>n = (2m + l)J2j + (2n + 1)Q 2 (35) 

now become complex when e 2 > 1 + c 2 . Indeed, in all of the examples studied in this 
section, the overall Hamiltonians are VT symmetric, and the transition to complex 
eigenvalues is a signal of the spontaneous breakdown of that symmetry. 
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4. Coupling to Generic Hermitian Hamiltonian 

In this section we examine the physical system described by a non-Hermitian VT- 
symmetric harmonic oscillator Hamiltonian Hi = p 2 + x 2 + 2ix coupled to a general 
Hermitian Hamiltonian H2 = p 2 + V(y). The only assumptions we will make are that 
H 2 is separately V and T symmetric. Thus, we assume that V(y) is real and is symmetric 
under parity reflection: V(y) = V(—y). Using only these assumptions we are able to 
show that the perturbative expansion for the ground-state energy is real up to 0(e 2 ). 
The ground-state eigenfunction 

rj( x ) = e -(*-M) a /2 ( 36 ) 

of the PT-symmetric harmonic oscillator satisfies the Schrodinger equation 

- rj"(x) + (x 2 + 2ix)r](x) = 2r](x), (37) 

whose ground-state energy is 2. Denoting the ground-state energy of H 2 by A, the 
ground-state wave function ip(y) satisfies the Schrodinger equation 

-f(y) + v(y)il>(v) = W{v). (38) 

We couple Hi and H 2 via the coupling term exy, so that the total Hamiltonian 
is H = Hi + H 2 + exy, with e considered as a small parameter. The ground-state 
eigenfunction $(x, y) of the combined system then satisfies the Schrodinger equation 

- 3> xx + {x 2 + 2ix)<S> - % y + V(y)$ + exy§ = E$. (39) 

Let us calculate E and y) perturbatively. The first three terms in the 
perturbation expansion for the energy are 

E = 2 + A + eE x + e 2 E 2 + ... (40) 

and we write 

$(x, y) = %(x, y) + e$i(x, y) + e 2 $ 2 (x, y) + . . . , (41) 

where 

M^y) = e~ ix+l)2/2 m- (42) 
The coefficient of e 1 in the expansion of Eq. fl39l) is 

- ($1)^ + {x 2 + 2ix)$i - ($!) ro + V(y)Q x = -xy® + E $i + E^ .(43) 

The solution to the homogeneous part of this equation is satisfied by $ (x,y) = 
T](x)ij}(y). Using the method of reduction of order, we therefore set $1 = $o( a; ) y)Q( x , y)- 
The integrating factor of the resulting equation is $o- Multiplying by this integrating 
factor gives the differential equation 

{®lQx)x + {®lQ y )y = {-xy + Ei)§l. (44) 

To find Ei we integrate this equation over all x and y and note that the integrals 
over the total derivatives vanish. This gives the following expression for Ei. 

f°° dxxe~( x+ ^ 2 f°° dyy<b 2 (y) 

771 J— 00 J— 00 a <* ~ ' (AK\ 

~ IZo dxe ~ {x+l)2 IZody^iy) ' 
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The integral over <p 2 (y) in the numerator vanishes because of parity symmetry. Thus 
Ei = 0. This result simplifies the differential equation satisfied by Q(x,y) to: 

(HQ*)* + (HQy)y = -xy*l (46) 
Proceeding to next order, we find the coefficient of e 2 in the expansion of Eq. (139]) : 
- ($ 2 ) M + {x 2 + 2zx)$ 2 - ($ 2 ) w + V{y)® 2 = -xy&x + E <t> 2 + £ 2 $ -(47) 

To solve this equation we again use reduction of order and set ^ 2 (x, y) = ®o(x, y)R(x, y). 
Multiplying by the integrating factor <E> , we get 

(§lR x ) x + ($ 2 Ry)y = (-xyQ + E 2 )$ l (48) 

The next correction to the energy comes from integrating this equation over all space: 

^ = I-oo dx I-00 d V xyQ(x,y)$&x,y) 

JZc dx fZc d yM x ,y) 2 

This time the correction does not vanish and we therefore must determine whether it is 
real or complex. Note that the integral in the denominator is real. 

Let us assume that the numerator / = dx dyxyQ(x, y)§l(x, y) is complex 
and expand &o(x, y) into its real and imaginary parts: 

$o(x, y) = ip 2 (y)e 1 ~ x2 [cos(2x) - % sin(2x)] = E(x, y) + iF(x, y). (50) 

Note that E is even in x and y, while F is odd in x and even in y. Next, we let 
Q = S + iT. Then flU reads 

[(E + iF)(S x + iT x )} x + [{E + tF)(S y + iT y )] y = xy(E + iF). (51) 
Taking the real and imaginary parts of this equation, we obtain 

(ES X - FT X ) X + (ES y - FT y )y = xyE (52) 
for the real part, and 

(ET X + FS X ) X + {ET y + FS y )y = xyF (53) 

for the imaginary part. We conclude that S is odd in x and y and that T is even in y 
and odd in x. Now 

dx / dyxy(ET + SF) . (54) 
-00 J —00 

Here ET is even in y, while FS is even in x. Thus, we have shown that I is real. 

This result establishes that for small enough e the ground-state energy remains 
real. However, we do not have sufficient information to determine whether there is a 
critical value of e at which the energy becomes complex. If, as in Ref. [19], we were 
able to establish that E is a Herglotz function of e, that is, that lm(E) has the same 
sign as Im(e), then we would indeed know that there was such a critical value because a 
Herglotz function that is entire must be linear [20], whereas we have shown that E 2 7^ 0. 
Unfortunately we are at the moment unable to construct a proof of the Herglotz property 
of E. 
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5. Summary 

We have shown in a number of examples that it is possible to couple Hermitian and 
'PT-symmetric non-Hermitian Hamiltonians together in such a way that the energy 
eigenvalues of the combined system remain real for sufficiently small values of the 
coupling e. In the matrix model and in all of the quadratic systems we have studied there 
is a critical range of the coupling, which, if exceeded, results in a complex spectrum. For 
coupling to a more generic PT-symmetric potential we have as yet no analytic proof of 
the existence of a critical point in e. 
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